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ABSTRACT 

We demonstrate that distinguishable gamma-ray burst pulses exhibit similar 
behaviors as evidenced by correlations among the observable pulse properties 
of duration, peak luminosity, fluence, spectral hardness, energy-dependent lag, 
and asymmetry. Long and Short burst pulses exhibit these behaviors, suggesting 
that a similar process is responsible for producing all GRB pulses. That these 
properties correlate in the observer's frame indicates that intrinsic correlations are 
strong enough to not be diluted into insignificance by the dispersion in distances 
and redshift. We show how all correlated pulse characteristics can be explained 
by hard-to-soft pulse evolution, and we demonstrate that "intensity tracking" 
pulses not having these properties are not single pulses; they instead appear 
to be composed of two or more overlapping hard-to-soft pulses. In order to 
better understand pulse characteristics, we recognize that hard-to-soft evolution 
provides a more accurate definition of a pulse than its intensity variation. This 
realization, coupled with the observation that pulses begin near-simultaneously 
across a wide range of energies, leads us to conclude that the observed pulse 
emission represents the energy decay resulting from an initial injection, and that 
one simple and as yet unspecified physical mechanism is likely to be responsible 
for all gamma-ray burst pulses regardless of the environment in which they form 
and, if GRBs originate from different progenitors, then of the progenitors that 
supply them with energy. 



Subject headings: gamma-ray bursts, methods: statistical and correlative studies 
of gamma-ray burst properties 
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Introduction 



Gamma-ray burst (GRB) prompt emission is generally representative of coherent, pulsed 
radiation, rather than of a random walk, a white noise process, or some other type of noise 
process. In recent years, pulse properties have provided increasingly valuable constraints on 
the physics responsible for GRB prompt emission. These properties have included (1) tempo- 
ral asymmetry characterized by longer decay t han rise rates, (2) hard-to-soft spectra l evolu - 
tion, and (3) broadening a t lower energies (e.g. Golerietskii et al. (Il983 ): Norris et al. ( 1996 ): 



Liang &: KargatisI Jl996[ ): iRamirez-Ruiz fc Fenimorel J200o[ l: iNorrisI J2OO2I ): iRvdel J2005l )l 



The vast majority of GRB pulses have similar shapes; at least they are sir nilar enough 



that most can be fitted by a simple four-parameter empirical pulse model (INorris et al. 
2OO5I ). T he model has allowed many isolated a n d semi-overlapping GRB pulses to be 



extracted flNorris et al.l l2005l : iHakkila et al.l I2OO8I : iHakkila fc Cumbed l2009al ). However, 
uniq ue extraction of heavily-ov erlapping pulses and very low signal-to-noise pulses is diffi- 
cult (IHakkila &: Cumbeell2009a[ ). and the characteristics of these unfitted pulses are unknown 
due to the dif ficulty in mak ing measurements as well as the inherent non-uniqueness of the 
solutions (e.g. lActonl (Il970l )). 



When applied to pulse data contained in individual energy channels, the pulse fit model 
demonstrates not only that pulses have l onger durations at l o wer energies, but that the y 
also peak later at these lower energies (IHakkila et al.l l2008l : IHakkila fc Cumbed l2009al ). 
This energy-dependent pulse lag is differeri t for each pulse (this can be verified using ei- 



ther the pulse peak lag (IHakkila et al 



20081 : IHakkila fc Cumbed l2009a] ) or by applying the 



cross- correlation function (iBandl Il997l ) to pulses, indicating that greater information con- 
tent resides in the pulse lags than can be obtained from measuring a single lag over the 
dur ation of the GRB. T he pulse peak lag inversely correlates with the pulse peak luminos- 
ity (IHakkila et al.l 120081 ). indicating that this relat ion is the basis of the lag vs. luminosity 
relation obtained for integrated prompt emission ( iNorris et al.ll2000l ). Furthermore, pulse 
duration also correlates directly with pulse peak lumin osity, allowing easily - measured pulse 
durations to be used as a GRB luminosity indicator ( Hakkila et al. 2008[). In addition to 



BATSE data, this t echnique has rec e ntly b een applied to HETE -2 ( lArimoto et al.l |2010[ ) 



and Swift data (e.g. IChincarini et al.l ( I2OIOI ): iMargutti et al.l (120101)'). and the technique ap- 



pears to provide a promising way of estimating GRB redshifts (IHakkila. Fragile, and Giblin 
2009bh . 



In recent preliminary work, IHakkila fc Cumbed ( l2009al ) have demonstrated that other 
pulse spectral and temporal characteristics (besides lag and duration) correlate with pulse 
luminosity, and thus with each other. Pulse duration correlates with pulse asymmetry and 
pulse fiuence and anti-correlates with pulse hardness and pulse peak fiux; all of these mea- 
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sures are thus pulse peak luminosity indicators. iHakkila &: Nemirofla ( 120091 ) have further 
shown that, within measurement uncertainties, pulse emission begins simultaneously at all 
energies. 

The ubiquitous presence of GRB pulses, along with their correlative properties, redirects 
attention from bulk GRB emission towards its constituent structures. Many luminosity- 
correlating properties attributed to bulk GRB pror npt emission are ac t ually p roperties of 
convolved pulses, and can be understood as such ( IHakkila fc Cumbed l2009a[ ). As previ- 
ously mentioned, a GRB's lag (INorris et al.l 120001 ) is a heterogeneous composite of the in- 



dividual pulse lags favoring the high intensity , short durati on pulses (jHakkila et al.l 12008 



Hakkila fc Cumbed l2009al : lArimoto et al.l |2010|) . Variability ( iReichart et al.l l200ll ) , measur- 
ing the changing luminosity of a GRB, increases and decreases with each pulse, and thus indi- 
cates a GRB's multi-pulse structure. Th e peak energy -Epea k (the peak of the vFy spectrum) 
correlates with the isotropic energy Eiso ( lAmati et al.ll2002l ). yet -Epeak is the time- integrated 
peak of the combined spectra of all pulses while E^so is the s ummed fluence luminosi ty of all 
pulses. The Internal Luminosity Function (ILF), or ip{L) ( iHorack fc Hakkilalll997l ). is the 
observed luminosity distribution within a GRB, and the quantity ip{L)AL represents the 
fraction of time that the GRB's luminosity is found between L and L + AL. The ILF sum- 
marizes information about the distribution of luminosity within a burst without explicitly 
describing the order in which the emission occurred, or if the emission was pulsed. However, 
the ILF measures a convolution of the luminos ity within each puls e along with the number of 
pulses which contribute to the total emission (IHakkila et al.ll2007| ). The correlation between 
GRB peak luminosities and ip{L) thus indicates that GRB luminosity is driven by pulse 
structure rather than being a generic property of the prompt emission. 

Since GRBs can be classified into Long and S hort bursts on the bas i s of duration. 



hardn ess, and fluence (e. g. iKouveliotou et al.l ( 19931): iMukherjee et al.l (Il998l ): IHakkila et al. 



( I2OO3I )). the preliminary IHakkila fc Cumbed (l2009al ) results demonstrate that the correla- 
tions between pulse properties are common to pulses originating in both the Short and Long 
GRB classes. Furthermore, there is substantial overlap between the properties of the shortest 
pulses in Long GRBs and Short GRB pulses, suggesting the existence of a single process. 

In this manuscript, we demonstrate that correlative GRB pulse properties are universal 
for Short and Long GRBs, and are thus progenitor-independent (if Short and Long GRBs 
originate from different progenitors). We further demonstrate that these simple pulse be- 
haviors place strong constraints on GRB physics. 
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2. Measuring Pulse Properties 



In this study, we analyze GRBs cataloged by the Burst And Transient Source Experi- 
ment (BATSE) on NASA's Compton Gamma-Ray Observatory (CGRO). By studying GRB 
pulses obtained from a single experiment, we limit our sample to a homogeneous one with 
well-understood instrumental and sampling biases. Inclusion of GRBs detected by other 
instruments (e.g. HETE-2, Swift, and Fermi) introduces a wide range of poorly-understood 
and pervasive systematic biases that complicate the study of GRBs across instruments: 
count rates from BATSE cannot be directly compared to those of other instruments be- 
cause each instrument is characterized by a different spectral range, temporal response, and 
energy-dependent sensitivity. For example, (1) Swift and HETE-2 are not sensitive at the 
higher energies of BATSE and Fermi, where the GRB photon count rates are highest, and 
(2) BATSE has a larger surface area and therefore a lower photon detection threshold than 
the similar detectors on Fermi's GBM experiment. BATSE's well-understood systematics 
have been extensively documented, which can be used to determine the viability of extracted 
pulse properties. 

The pulses described here are part of a BATSE GRB pulse catalog currently being 
produced (Hakkila et al., in preparation) by an analysis of the Current BATSE Cata- 
log (|http : //www . batse . msf c . nasa . gov/batse/grb/catalog/ current/ ) . The process in- 
volves extracti ng pulses from summed 4-channel 64 ms GRB data using the Bayesian Blocks 
methodology (IScarglelll998l ) and subsequently fi tting the possible pulse time intervals with 
the four-parameter model of iNorris et al.l ( 120051 ). 



The Bayesian Blocks process (jScarglelll998l ) at the root of this analysis is a statistical 
model by which a photon-counting light curve is segmented into a sequence of discontinuous 
constant levels; each level indicates that the background is statistically different than the 
segments preceding and following it. The result of a Bayesian Block segmentation looks like 
histogram with variable-width bins. Segments with elevated backgrounds are regions that 
possibly contain pulses. 

We assume that cand idate pulses can be modeled with the empirical functional pulse 
form of iNorris et al.l (|2005[ ): 

I{t) = AAexp[-"i/(*-*=)-(*-*^)/"^l, (1) 

where t is time since trigger, A is the pulse amplitude, tg is the pulse start time, ti and T2 
are characteristics of the pulse rise and pulse decay, and the constant A = exp [2(ri/r2)^/^]. 
The 4-channel pulse peak time occurs at time Tpeak = tg + a/t^. As each segment is 
searched for a pulse, statistically-insignificant intervals are merged into the surrounding 
intervals until an optimal set of model pulses is identified for the GRB. The criterion for 
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i dentifying a statist ically-significant pulse is based on the idea of a dual timescale peak flux 



( iHakkila et al.ll2003l ). which does not bias long, faint pulses relative to short, bright ones. The 
long timescale peak fiux is based on the 4-channel fiuence (time-integrated photon counting 
fiux S'tot) divided by the pulse duration w while the short timescale is based on the 4-channel 
64 ms peak fiux (^64)- 

The criterion generally used in our analysis to select pulses is 

P64Stot/w>n% (2) 

where b is the average background counts in a 64-ms interval, w is the pulse duration 
(described below), and n is the number of standard deviations above background that the 
pulse peak needs to be in order to be considered. The default setting of the automated code 
is n = 3, which works well for most isolated pulses and many partially overlapping pulses. 
However, operator assistance is needed when inspection indicates that distinct pulses have 
not been properly fitted. For example, thresholds as low as n = 1 have been required to 
fit some long, faint, isolated pulses for which the BATSE background rate is particularly 
well-behaved. Larger values of n (e.g. n > 10) are needed to extract some overlapping 
pulses in complex GRBs because the greatest source of noise to fitting a pulse can be the 
local background caused by overlapping pulses. 

Four examples of this process are demonstrated in Figure 1; single-pulsed BATSE trig- 
gers 1221 and 432, and multi-pulsed triggers 1443 and 1468. Each Bayesian block segment is 
identified by vertical lines which overlay the raw counts data. The pulses fitted to the data 
are also shown, as is the overall fit. 

Pulses fitted in the summed 4-channel data are subsequently extracted from individual 
energy channels using the 4-channel pulse parameters as starting points. Although the 
process is fairly robust, low intensity pulses, pulses that strongly overlap, and/or very short 
pulses indistinguishable from Poisson noise are not considered to have been successfully 
fitted, and are excluded from this discussion. 

A variety of model-dependent, observable pulse properties are extracted from the mul- 
tichannel parameters tg, A, ti, and T2. These observables are extracted independently for 
energy channels 1, 2, 3, and 4 as well as for the summed BATSE 4-channel data, although 
low signal-to-noise makes channel 4 pulse properties very difficult to measure. These observ- 
ables (chosen so that they can be compared to properties of the prompt emission) include 
pulse durations, pulse peak lags, pulse peak fiuxes, pulse asymmetries, and pulse spectral 
hardnesses. 



The pulse duration w is obtained from the summed four-channel data and is defined 
as w = [9 + Vl^ txIt'^I'^; this is the interval between times when the pulse amplitude is 
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Ae~^ (note that iNorris et al.l (120051 ) uses the shorter interval between times when the pulse 
amplitude is Ae~^). Pulse peak lags I are the differences between the pulse peak times in 
different energy channels (pulse peak times are given by Tpeak = ts + a/tTt^)- Pulse peak lags 
can be obtained for any pulse between two energy channels, although we define the standard 
pulse peak lag / as that measured between energies of 100 to 300 keV (BATSE channel 3) 
and 25 to 50 keV (BATSE channel 1). Pulse durations and pulse peak lags are measured in 
units of seconds. 

The pulse peak flux p2m is the peak flux (in counts cm"^s~^) of the fitted pulse defined 
on the 256 ms timescale. The pulse fluence (S) is the time-integrated flux summed over the 
pulse duration in the 25 to 300 keV range (channels 1 through 3), given in units of ergs cm~^. 
We define the fluence in this way (rather than in counts cm~^) so that the results can be more 
directly compared to published BATSE burst fluences. In order to make this conversion, we 
have assumed an approximate average energy for each BATSE channel, and we have also 
excluded channel 4 fluence because of BATSE's lower quality instrumental response at these 
energies. 

The pulse asymmetry k is defined as k = w/(3 + I^txIt-i); it ranges from a value of 
K = for a symmetric pulse to k = 1 for an asymmetric pulse with a rapid rise and slow 
decay. Hardness ratios HR are constructed by dividing pulse fluences in two different energy 
channels; we use the hardness ratio HR= S^/Si, where is the channel 3 energy fluence 
and 5*1 is the channel 1 energy fluence. Both the pulse asymmetry and the hardness ratio 
are unitless parameters. 

The GRB s have been classified using the BATSE-dependent data mining definition of 



Hakkila et al.l ( 120031 ) . which classifies GRBs as Short if they satisfy the inequality (T90 



< 1.954) OR (1.954 < T90 < 4.672 AND HR321 > 3.01), where HR321 = ^3/(^2 + ^i) 



(IMukherjee et al.l 119981 ). and where 5*2 is the channel 2 fluence. If GRBs do not satisfy the 



aforementioned inequality, then they are classified as Long. 

Electronic Table 1 contains the GRB pulse properties used in this study, along with 
their uncertainties. This table contains preliminary data from the BATSE pulse catalog 
currently under development. Current pulse extraction has yielded to date 1338 pulses in 
610 BATSE GRBs. 



3. Pulse Property Correlations 



The properties of the GRB pulses examined in this study exhibit many strong correla- 
tions; these can be seen in mutipanel Figure 1. This figure demonstrates the relationships 
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between the observable pulse properties of duration, lag, peak flux, hardness ratio, asym- 
metry, and fluence. Pulses from the Short class of GRBs are shown in blue, while pulses 
from the Long GRB class are denoted in red. Many correlations and anti-correlations are 
apparent by inspection, and the statistical significances of the correlations are summarized in 
Table 2. Since Short GRB pulses have short durations and large amplitudes, they occupy an 
extreme position in each distribut ion. The corre l ations and anti-correlations illustrated here 
have been described previously in lHakkila et al.l ( l2008l ) and lHakkila &: Cumbed ( l2009al ). but 
are much more statistically significant given the larger size of this pulse sample. 

A strong correlation exists between pulse duration and pulse peak lag; this correlation 
is present for both Long and Short burst pulses. However, not all pulse lags are positive; 
one-third are zero or negative. Most of the negative lags are consistent with their mea- 
surement uncertainties (indicating that they are simply too small to accurately measure), 
although large negative lag values are suggestive of systematic ex traction issues involving 
unresolved overlapping pulses (e.g. see lHakkila fc Cumbed (l2009al )). Figure 2 demonstrates 
the relationship between pulse duration and pulse peak lag for the pulses with positive lags. 
Assuming that pulses have positive lags, pulse peak lag appears to be a proxy for hard-to-soft 
pulse evolution. By this we mean that the timescale for pulse emission correlates with the 
timescale of the hard-to-soft evolution. This effect is not distorted by cosmological redshift z 
to first order, because the values of both duration and lag in the emitted frame are (1 + z)~^ 
the values measured in the observer's frame (in making this statement, we purposefully ig- 
nore effects caused by the fc-correction and the spectral range of the instrument, which will 
be discussed later). 

Pulse duration anti-correlates with pulse peak flux for both Long and Short burst pulses. 
The brightest pulses (as measured by peak flux) found within observed GRBs are also the 
shortest pulses. This is surprising given the strong c o rrelation betwe e n pul s e peak luminosity 



and p ulse duration (e.g. iHakkila fc Cumbed (l2009al ): lArimoto et al.l (120101 ) : IChincarini et al. 



(I2OIOI )). because it indicates the correlations are strong enough to not be diluted into in- 
significance by the dispersion in GRB distances. The photon peak flux on the proper motion 
distance du and redshift z is proportional to [(1 + 2;)(i|^]~^. Again, the first-order effects 
of time dilation cancel in the comparison between duration and peak flux, so that only the 
inverse square law disperses the relation. However, the similarities between the peak flux 
vs. duration distribution and the luminosity vs. duration distribution suggest that the peak 
flux distribution of GRB pulses is closely related to the pulse peak luminosity function. We 
note that the pulses from Short GRBs appear to extend the correlation found for Long burst 
pulses such that the shortest pulses from Long GRBs overlap the Short GRB distribution. 



Pulse spectral hardness is anti-correlated with pulse duration such that shorter, brighter 
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pulses are generally harder than longer, fainter pulses. The relationship is weaker for Short 
GRB pulses than for Long ones, although it still seems significant. We note that spectral 
hardness has no z-dependence, but it does require a fc-correction, which may be significant 
for bursts that differ in redshift by factors of order 2. Since duration has not z-dependence, 
there is a scatter in the correlation caused by the redshift dependence of the duration relative 
to the spectral hardness. 

Pulse asymmetry also correlates with pulse duration, such that the most asymmetric 
pulses are long duration, long lag, low peak flux, soft pulses. Since pulse shape is difficult to 
measure for pulses having durations approaching the 64 ms timescale (Poisson background 
fluctuations can significantly mask pulse characteristics). Short GRB pulses tend to be fitted 
with a bimodal asymmetry distribution with asymmetries approaching the extreme values 
of either k = or k = 1, and show no statistically significant correlation with pulse dura- 
tion. Pulse asymmetry also has no z-dependence, so correlations between it and other pulse 
properties is dispersed in those properties by a factor of (1 + z)^^. 

Pulse energy fiuence S is strongly correlated with pulse duration; longer pulses contain 
more energy than short pulses. A strong anti-correlation also exists between pulse peak 
flux and pulse fiuence, indicating that the pulses that emit their energies in the shortest 
time interval generally produce the fewest photons over their duration. And, to a first 
order approximation, the harder spectrum of a short pulse does not provide enough energy 
to compensate for the larger number of photons produced by a longer pulse. Thus, the 
pulses with the largest energy fiuence are long, soft, low peak flux pulses. Because we have 
chosen to measure fiuence i n energy units, t he red shift dependence is again proportional to 



[(1 + z)d1,j] ^ (e.g. see also iMeszaros et al.l (120111 )). Thus, redshift only causes a dispersion 



in S when compared to k, and HR, which have no redshift dependence. However, like peak 
flux, the inverse square law introduces a d~j^ dispersion when being compared to all pulse 
properties other than peak flux. 

Most of the remaining relationships between GRB pulse properties described in Table 2 
are as expected, given the aforementioned correlations and anti-correlations (the exception 
is the hardness vs. fiuence relation, which will be discussed in section 5). The difficulty in 
accurately measuring lag and asymmetry for very short pulses makes it hard to accurately 
measure correlations between these properties and other ones for short duration GRB pulses. 

We note that, since lags and durations of Long GRB pulses are indicators of pulse peak 
luminosity, the observer frame attributes of spectral hardness, pulse asymmetry, fiuence, 
and peak flux are also luminosity indicators. We also note that pulses from Long and Short 
GRBs appear to represent a continuum of values; we will return to this point in section 8. 



-9- 



4. Intra-Burst Pulse Properties 

The previous discussion refers to correlations between all pulses, regardless of the bursts 
in which they originate (e.g. inter-burst pulse correlations). We can also look at the distri- 
butional properties of pulses within GRBs (e.g. intra-burst pulse properties). Although we 
have only had limited success fitting pulses in complex GRBs, we have many bursts with 
two, three, and four pulses in them. However, since only a small number of Short GRBs 
have multiple distinct pulses on the 64-ms timescale, we treat multi-pulsed Short GRBs as 
if they are part of the same sample as the multi-pulsed Long GRBs. 

Averaging together the properties of pulses within multi-pulsed GRBs gives us some 
sense of the bulk intra-burst pulse distribution. Figures 4 through 7 indicate the mean values 
of several key, easily-measured pulse properties measured within individual GRBs: duration, 
peak flux, hardness, and fluence (the error bars indicate the mean error for each averaged 
measurement). The m,ean pulse properties for multi-pulsed GRBs correlate similarly to the 
bulk properties of individual pulses. In other words, average pulse duration is apparently 
a good indicator of average pulse brightness, hardness, lag, and fluence, just as is found 
in single-pulsed bursts (although the variation of pulse asymmetries, coupled with large 
measurement uncertainties, make them an uninformative average pulse property). 

It is more difficult to ascertain if the pulses within multi-pulsed GRBs themselves obey 
the correlations identified in Figure 2 and Table 2. Overlapping pulses in multi-pulsed 
bursts are the most difficult to fit and typically provide the least reliable pulse property 
measurements. In other words, bursts with the largest numbers of pulses tend to have the 
smallest numbers of accurately fitted pulses. Thus, in attempting to verify the correlative 
relationships between pulses within bursts, we will be limited to intra-pulse data consisting 
of bursts with a small number of pulses having well-measured properties along with bursts 
having many pulses but with less reliable pulse property measurements. 

As a simple test of whether the correlations exist, we calculate the linear gradients 
(slopes) of pulse parameters within multi-pulse GRBs. We apply least squares to the prop- 
erties of all pulses fitted within each multi-pulsed GRB to learn more about the intra-burst 
pulse distribution. The resultant gradient of each pulse property relative to the others in- 
dicates whether the properties within the burst are correlated (positive) or anti-correlated 
(negative). The percentage of multi-pulsed GRBs exhibiting correlated behaviors is then 
compared to the results of the inter-burst pulse duration shown in Figure 2 and Table 2 
and is demonstrated in Table 3; if somewhat more than 50% of the multi-pulsed GRBs have 
the same general gradient as that obtained from the inter-burst pulse distribution then the 
two distributions have the same general behavior. Since these percentages average 62%, the 
behaviors of the intra-burst pulse distribution have the same sense but otherwise appear to 
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be generally weaker than the power-law forms of the inter-burst pulse distribution. Although 
these weaker relationships might indicate weaker intrinsic relationships between pulse prop- 
erties within a burst, they might simple reflect a greater difficulty in accurately measuring 
these properties in complex GRBs. 



5. Possible Contributions of Sampling Biases to Pulse Property Correlations 

The observed pulse correl a.tions are not caused by the wa y pulses have been ext racted, or 



by other selection biases (e.g. iHigdon fc Lingenfelterl (jl996[ ): iMeegan et al.l (|2000[ )). Figure 



8 shows that the peak flux vs. duration correlation is strongly present even for single-pulsed 
GRBs triggering on the 256 ms timescale (from the BATSE Cmax/C'min Table found at 
|http : / / www ■ batse . msf c . nasa . gov/batse/grb/ catalog/ curre nt/ [ ) . All of the pulses in- 
cluded in Figure 8 have been selected on the basis of peak flux rather than the dual-timescale 
pulse trigger used in the selection of pulses from multi-pulsed bursts. The sample is lacking 
in short pulses that trigger preferentially on the 64 ms timescale (and at the cost of long 
duration, low peak flux pulses) and in long pulses that trigger preferentially on the 1024 
ms timescale (at the cost of short duration pulses), yet still shows a pronounced correlation 
between peak flux and duration. The overall sample shown in Figure 2 has an additional 
advantage of triggering on BATSE's three different trigger timescales, since this sample in- 
cludes pulses that triggered on 64 ms OR 256 ms OR 1024 ms. The multi-timescale trigger of 
the entire sample is thus less susceptible to selection biases than a catalog selected solely on 
the basis of a single timescale trigger would be. As the remaining pulses in the pulse catalog 
have also been selected using a dual-timescale trigger (peak flux vs. fiuence), it appears that 
biases related to the trigger are not responsible for the existence of the pulse peak flux vs. 
duration correlation. 

We can also determine if the dual timescale peak flux used to select pulses is some- 
how responsible for creating the observed pulse property correlations. We have estimated 
the threshold dual timescale peak flux limits of extracted GRB pulses using equation 2 in 
conjunction with mean BATSE background counts; the results are shown in Figure 9 in 
terms of the two constituent components of the dual timescale peak flux: the counts fiuence 
divided by the duration, and the 64 ms peak flux. The solid line represents the standard 
dual timescale peak flux threshold described in equation 2 when 5 = 3, while the dotted 
line represents the near minimum threshold when 6 = 1. The approximate durations of the 
plotted pulses are indicated on the diagram, since pulse duration correlates strongly with 
fiuence. The diagram suggests that background is not an issue for pulses with durations 
shorter than ^ 10 s, although it makes it increasingly difficult to recognize and extract 
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pulses longer than 100 s. Thus, because the slope of the data is roughly perpendicular to 
the pulse-fitting threshold, the threshold is not responsible for the strong pulse property cor- 
relations observed in short, bright pulses. Furthermore, correlations continue to adequately 
describe peak fluxes and fiuences of pulses fainter than this threshhold, suggesting that the 
correlations also hold for longer, fainter pulses. Although we have demonstrated that trigger 
biases are not responsible for the pulse correlations, it does appear that they prevent us 
from detecting very long, low-amplitude, single-pulsed bursts (e.g. pulses on the order of 
hundreds to thousands of seconds) . 

The analysis has some other important selection biases as well. As mentioned before, 
the analysis purposefully excludes bright pulses which cannot be uniquely identified due to 
pulse overlap. Similarly, pulses that are too faint and/or too short (if there are any of these) 
will be excluded on the basis that they can be confused with Poisson noise. However, the 
evidence indicates there does not seem to be a systematic exclusion of pulses capable of 
explaining the observed pulse correlations. 

The one statistically-significant correlation that represents a difference between Long 
and Short GRB pulses is that between hardness and fluence (see Table 2). Short GRB 
pulses exhibit a strong correlation between hardness and fluence, whereas Long GRB pulses 
exhibit a moderately strong anti-correlation between these properties. This disparity is not 
real, but is due the fluence deflnition being used in this analysis. Short duration pulses, 
and particularly those belonging to the Short GRB class, arc very hard and often have 
a signiflcant amount of channel 4 fluence. Since our fluence deflnition excludes channel 
4 (because most pulses are too difficult to accurately fit at these energies), and since the 
average energy per photon is largest in channel 4, we are systematically underestimating the 
energy contribution of the shortest, hardest pulses. The shorter and brighter a pulse is, the 
more our fluence definition has underestimated the high energy contribution. 

As has been mentioned, the effect of redshift cancels in almost all the aforementioned 
pulse property correlations to first order. However, energy-dependent corrections resulting 
from the instrumental bandwidth (/c-corrections) are also important, although more difficult 
to assess. Pulses observed at large distances should have their high energy emission shifted 
into low-energy channels and their low-energy emission shifted out of the instrumental pass- 
band. Increased distance will thus cause -Epeak to shift to lower energy channels, sometimes 
moving it out of the instrumental passband. Most distant pulses will appear softer as -Epcak 
is shifted to lower energies. Since pulses have shorter durations at high energies than at low 
energies, highly-redshifted pulses should appear to be shorter than their nearby counterparts; 
this effect works in the opposite sense of time dilation. Instruments that are insensitive to 
high energy emission, such as Swift and HETE-2, cannot generally detect the high energy 
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pulse peaks occurring early in the pulse emission; pulses from these experiments thus often 
appear to be less peaked and more time-symmetric than their BATSE counterparts. 



Pulse Evolution 



As evidenced in Table 2, the correlated characteristics of GRB pulse emission are all 
related to the pulse spectral lag, which measures hard-to-soft pulse evolution. The pulse peak 
lag, denoting the pulse peak delay between BATSE channel 3 and BATSE channel 1, is only 
one measure of this internal pulse evolution. Pulse peak lags also exist between all BATSE 
energy channels, reflecting the fact that the pulse intensity peaks first in the highest energy 
channel and subsequently in progressively lower en ergy channels. Furthermo re, pulses start 
near-simultaneously in all BATSE energy channels (IHakkila fc Nemirofj|2009t ): when coupled 
with the fact that pulses peak earlier at high energies, this indicates that hard emission 
dominates over soft emission ea rly in the pulse, and d ecays away faster to be replaced by 
soft emission later in the pulse (IRichardson et al.l 119961 ). This evolution is demonstrated in 
Figure 10 for single-pulsed BATSE Trigger 1883. 

There are other techniques for identifying pulse evolution. Perhaps the mos t direct way 



is via the time variability of the uF^, pulse spectrum via the Band spectral model (IBand et al. 
19931 ). Three characteristics are used to parameterize the spectral shape: the low-energy 



power-law index a, the high-energy power-law index /3, and the peak of the spectral energy 
distribution -Epeak- Pulse spectral softening indicates a systematic evolution of these three 
spectral parameters, and primarily a shift from high -Epeak to low -Epeak- 

Spectra of hard-to-soft (classic FRED, denoting Fast Rise Exponential E >ecay) pulses 



have b een observed to unde r go a decay of . Epeak t hroughout t he pul se (e.g. iCrider et al. 
fll999h : lKocevski et al.l J2003h : iKaneko et al.l J2006h : IPeng et all J2009ah l The hardest spec- 
tra are observed early in pulse evolution, with the highest -Epcak values occurring when pulse 
flux is just beginning to rise. Pulses undergo a time decay of -Epeak that appears to b e 
exponential, but which can be fit by a linear function in many cases (ICrider et al.lll999l ). 
Interestingly, there is nothing special about the -Epeak at the time of the pulse peak flux; 
there is no change in the spectral decay that delineates the pulse peak from any other time 
during pulse spectral evolution. 



Although hard-to-soft evolution is the most common pulse evolution (e.g. ICrider et al. 



(Il999[). intensity track i ng pulse evolution has be e n identified in a sma l ler number of pulse s 



Ford et al.l (119951 ): iLiang fc Kargatid (119961 ): iKaneko et all (120061 ): IPeng et al.l fl2009bh : 



Lu et all torn ). In intensity tracking evolution, spectral hardness is found to correlate with 
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pulse flux. 

Intensity tracking pulse evolution appears to be inconsistent with pulse spectral lag 
measurements, since it cannot result in lags that extend across multiple energy channels. 
To better understand intensity tracking pulses, we select a sample of 22 hard-to-soft and 



Pene: et al. 


f 2009a 


b); 


Lu et al. 


(2010) 



In Table 3 we summarize the characteristics of these pulses and the GRBs from which 
they are extracted. Six pulses must be excluded from our analysis, as they have been classified 
as both intensity tracking and hard-to-soft in various studies, leaving eleven hard-to-soft 
and five intensity tracking pulses. Three of the intensity tracking and two of the hard-to- 
soft pulses are discarded because they have been fitted in our pulse catalog by two or more 
overlapping pulses. Moreover, all but one of the pulses classified as intensity tracking are 
found to originate in complex, multi-pulsed GRBs, where overlapping pulses are common. 
Hard-to-soft pulses are more often found as isolated pulses. We discard these ambiguous 
pulses from our sample, reducing it to nine hard-to-soft and two intensity tracking pulses 
that satisfy our criterion in equation 2. Four of the pulses (intensity tracking pulse 1733 and 
hard-to-soft pulses 2662 and 2665) have bumps prior to and/or after the pulse peak strongly 
suggesting the presence of additional long duration, low amplitude pulses. The remaining 
intensity tracking pulse (the first pulse in complex BATSE trigger 2812) appears to have two 
peaks that are most readily apparent in channels 2 and 3. Reconsidering these four pulses on 
the basis that they are potentially ambiguous, our sample contains only seven hard-to-soft 
pulses and no intensity tracking pulses. 

Although the evidence preliminarily suggests that intensity tracking pulses are composed 
of overlapping hard-to-soft pulses, we re-examine some of our excluded pulses to see if there 
is some characteristic that they have in common. Indeed, one of the intensity tracking pulses 
and two of the pulses ambiguously classified have similar long duration, time-symmetric light 
curves. These pulses all have large non-Poisson scatter occurring near-simultaneously in sev- 
eral energy channels, which is consistent with the characteristics defining intensity tracking 
pulses. However, the instantaneous -Epcak of these pulses starts high at the pulse onset before 
it begins to track the flux, and -Epeak resets each time that non-Poisson scatter occurs. In 
other words, Occam's razor suggests that each of these pulses is actually composed of over- 
lapping, short duration, hard pulses that are themselves consistent with the aforementioned 
correlated pulse properties. Thus, intensity tracking pulse behavior is unlikely, and simple 
hard-to-soft pulse behavior appears to be the dominant and perhaps only mode of pulse 
evolution. 



The presence of faint, unresolved pulses provides a possible explanation for many "out- 
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liers" shown in Figures 2 and 3; despite our precautions, pulses exist exhibiting large disper- 
sions from the otherwise strong correlative relationships. We believe that these are likely to 
be unresolved, merged pulses. Preliminary Monte Carlo analysis suggests that it is difficult 
to resolve two similar pulses peaking near-simultaneously, and that the characteristics taken 
on by the merged pulse are not simply a sum of the constituent pulse characteristics. 



Generally, these results suggest that the empirical four-parameter pulse model of lNorris et al. 



( 120051 ) can adequately fit the vast majority of resolved, non-overapping pulses. 



7. Peak Luminosity: a Consequence of Pulse Evolution 

Pulse duration and lag strongly anti-correlate with pulse peak luminosity. They also 
anti-correlate with pulse peak flux. The relationship between pulse peak flux and pulse 
duration is similar to, but has a broader spread than the relationsh i p bet ween pulse peak 



luminosity and pulse duration (compare Figure 2 with iHakkila et al.l ( l2008l ) Figure 4). This 
is somewhat surprising, since the pulse duration vs. peak luminosity relation was established 
using GRBs with known redshift and corrected for luminosity distance, the inverse square 
law, and cosmological time dilation. That ant i- correlations between pulse peak flux and 
duration are apparent without making these corrections (and without even converting photon 
flux to energy flux) indicates that dispersive effects caused by the inverse square law and 
relativistic cosmology do not cause enough observational scatter to destroy the intrinsic 
relationships. This is contrary to common understanding, as relativistic cosmology was once 
thought to be of central i mpor tance in the understanding of observed GRB characteristics 



(see the review by iPiranI (120051 )). 



It is perhaps not surprising that intrinsic pulse characteristics are so much stronger 
than relativistic cosmological effects given the presumed relativistic ejecta origins of GRB 
prompt emission. GRB pulse durations span flve orders of magnitude in the observer's 
frame, yet time dilation should stretch pulse durations by at most a factor of ten (for z = 9). 
Presumably, other intrinsic effects also dominate. For example, relativistic effects caused 
by the Lorentz factor P within the GRB jet {e.g. enhanced brightness due to relativistic 
beaming) can dominate because the blueshift of the jetted material (with P between 100 
and 1000) greatly exceeds the cosmological redshift. These effects could explain why bright 
GRBs tend to be dominated by complex light curves containing multiple, short- duration, 
short-lag, luminous pulses, while the faintest GRBs more often have long duration, long lag 
single pulses. Because pulse duration also correlates with fluence, hardness, and anisotropy, 
each of these observed pulse characteristics are also pulse luminosity indicators. 
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Our evidence suggests that the pulse peak flux is merely one step in the Epeak decay 
of a pulse spectrum; it does not appear to indicate a phase in the pulse spectral evolution 
that is different from any other. This is non-intuitive, because peak flux is a pinnacle in 
the light output, and therefore seems to be a transition between the pulse rise and the pulse 
decay. However, Figure 10 demonstrates that the pulse peak shifts from high to low energies, 
and as such is susceptible to the spectral response of the instrument (the instrument cannot 
measure the early contribution to the pulse peak at energies higher than its upper energy 
limit, nor can it measure the late contribution at energies lower than its lower energy limit). 
Despite this limitation, peak flux correlates with all other observable pulse properties. How 
does the detector sample observable pulse properties so that they always correlate with pulse 
peak flux, since the value of the pulse peak flux depends on the instrumental response as 
well as the redshift of the GRB and the intrinsic hardness of the pulse? 

This problem is more easily explained by using Occam's Razor to recognize that common 
(and perhaps ubiquitous) hard-to-soft evolution defines a pulse better than the increase and 
subsequent decrease in flux. Once we recognize this, the peak flux is seen just to be a 
byproduct of the pulse's spectral evolution, rather than its climax. The peak flux therefore 
represents the predictable sum of flux contributions from different energy ranges produced 
by the pulse energy decay. By understanding the mechanism capable of producing the pulse 
spectral evolution, we will understand why all bulk properties of the pulse are correlated. 



8. The Universality of Long and Short GRB Pulses 



The data presented in Figures 2 and 3 and summarized in Table 2 suggests that the 
observable properties of GRB pulses, regardless of how the bursts have been classified, make 
up a single distribution. The observable properties of Long and Short GRBs correlate in 
similar ways and overlap such that the shortest pulses of Long GRBs have characteristics 
indistinguishable from typical pulses in Short GRBs. In other words, there seems to be 
a mechanism capable of producing GRB pulses with correlated properties regardless of the 
environment in which they form and potentially of the progenitors from which they originate. 

The pulse peak flux vs. duration relation supports the idea that Long and Short GRB 
pulses constitute a single distribution, independent of pro genitor or environm ent. A recent 
analysis of Short GRBs observed by the Swift exper iment (iNorris et al.ll201ll ) demonstrates 



Norris et al 



(1201 If ) Table 1), and the 



that this relation continues to ms timescales (see 
relation holds regardless of whether the GRB containing them exhibits extended emission 
(EE) or does not exhibit extended emission (non-EE). 
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The overlap between Short and Long GRB pulse properties is substantially larger than 
indicated in Figures 2 and 3 because few of the shortest duration Long GRB pulses have 
been included in this analysis. These pulses are often found in complex GRBs overlapping 
other pulses, and thus cannot be fitted easily or uniquely. 

Is it possible on the basis of luminosity that Short and Long GRB pulses do not consti- 
tute a single distribution? Although other correlated GRB properties appear to be part of a 
continuous distribution, Short GRBs are known to be less luminous than Long GRBs, sug- 
gesting that the peak flux vs. duration relation is the result of a serendipitous confluence of 
conditions such as jet beaming angles, jet Lorentz factors, and burst redshifts that give Short 
and Long pulses similar peak fluxes when in reality their luminosities are quite disparate. 
To test this hypothesis, we take a sample containing po st-BATSE Short and Long GRBs 



with known isotropic luminosities ( IGhirlanda et al.ll2009l ). Although we have not explicitly 



measured the properties of these pulses using our pulse fitting method, we can still roughly 
estimate the average isotropic luminosity per pulse by dividing the GRB luminosity by the 
number of pulse peaks. The results, shown in Figure 11, indicate that Short and Long GRB 
pulses indeed have comparable luminosities, even though they originate from bursts with 
different luminosities. Long GRBs appear to be more luminous than Short GRBs because 
they are composed of more pulses having generally larger energy fluences. 

It is difficult to accurately measure pulse properties other than peak intensity and du- 
ration for Short GRBs due to the short timescale on which these properties are measured. 
However, pulse peak intensities, durations, and luminosities provide evidence beyond the 
pulse property correlations suggesting that Short and Long GRB pulses constitute a single 
phenomenon. If Short and Long GRBs indeed originate from different progenitors, then 
these pulses suggest that they are a progenitor-independent and environment-independent 
phenomenon. 



9. Conclusions 

Pulses are the basic building blocks of GRB prompt emission. However, pulse properties 
are challenging to measure: they exhibit large individual variations due to intrinsic differences 
(e.g. intrinsic -Epeak), systematic issues with background rates and instrumental response 
characteristics, and pervasive pulse overlap. Despite this challenge, all observable GRB pulse 
properties appear to correlate statistically, with correlated pulse properties being strong in 
the observer's frame. This surprising result indicates that the inverse-square law and some 
effects of relativistic cosmology are only of secondary importance to the internal physics 
responsible for producing the intrinsic pulse properties. GRB bulk properties are constructed 
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by combining and smearing out pulse characteristics in ways that potentially lose valuable 
information. Since pulse observables have been measured from a four parameter model, it 
is not surprising that the pulse model parameters themselves are highly correlated. The 
correlated nature of the observables constrains the underlying physical process capable of 
producing pulses: the process cannot have too many free parameters because the observed 
pulse properties might then be largely uncorrelated. Pulses start near-simultaneously in 
all energy channels and undergo spectral softening after this initial injection. Furthermore, 
it appears that pulses exhibit remarkably similar temporal evolutions as is indicated by 
their similar shapes. GRB pulses have progenitor-independent correlated properties: pulse 
physics appears ubiquitous and is easily replicated across a variety of GRB environments. 
Since Long GRBs can have high flux, short duration pulses, GRB pulse properties are a 
red herring in GRB classification; the internal pulse duration distribution may be a more 
valuable classification tool. The application of new statistical techniques (T. Loredo, private 
communication) will allow for better understanding of pulse-measuring limitations in GRBs 
characterized by considerable pulse overlap. 

The correlative properties and ubiquitous nature of GRB pulses places strong constraints 
on GRB physics. The standard model of GRB prompt emission involves kinematic energy 
injection via relativistic jets, after whic h the medium cools. T he standard spectral model 
produces either a synchrotron spe c trum ( 
spectrum (e 



Goodman (1986) 



Rees &: Meszaroslll994). or a thermal plus po A ver law 



Meszaros fc ReesI J200oh : lOaigne fc MochkovitchI j^OoJ; 



Rydd ( 12 004 |2005[ ): iGuiriec et al.l ( 120111 )). These models, and all other viable models, need 
to explai n the corre l ated p ulse properties and pulse evolution described in this paper. For 



example, iBoci et al.l ( 120101 ) point out that the correlative pulse relations cannot be obtained 
as a direct and simple consequence of the standard synchrotron shock model, since this has 
no time-dependent component. Many temporal and spectral pulse properties are in rough 
agreement with predictions related to emitting region curvature within the relativistic jet, 
although these are mo del- dependent and rely on offset viewing angl e to the jet, time-evolutio n 
of the physical process, and other model-dep e ndent effects ('e.g . Iloka fc Nakamural (1200 11) : 
Lu et al.l (l2006l ): lHafizi fc MochkovitchI (I2OO7I ): iBogi et al.l ( |201ol ): krimoto et al.l (I2OI0I ) ~ 
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Note. — Electronic Table 1 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. Table 1 contains the pulse properties for 1338 pulses from 610 BATSE GRBs used 
in this study, along with their uncertainties. These are preliminary data from the BATSE pulse catalog currently under development. 
Column 1 contains the pulse number (BATSE trigger ID plus fitted pulse number) , column 2 contains the burst classification (S=Short 
and L=Long), columns 3 and 4 contain the duration w and duration uncertainty (Tu, (measured in s), columns 5 and 6 contain the 
lag I and lag uncertainty cr; (s), columns 7 and 8 contain the 64 ms peak flux P64 E^nd peak flux uncertainty (Tp64 (counts cm^'^s"-'-), 
columns 9 and 10 contain the energy hardness ratio HR and HR uncertainty ithRi column 11 contains the asymmetry k (uncertainties 
(Tk all assumed to be ±0.05), and columns 12 and 13 contain the channel 1 -h 2 H- 3 energy fiuence S and fiuence uncertainty as (ergs 
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Pulses from Long GRBs 
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Pulses from Short GRBs 
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Table 2: GRB pulse correlations. Spearman rank order correlation probabilities that the 
pulse characteristics in question are as expected from random, uncorrected distributions. 
Anti-correlations are indicated by a f, and boldface numbers indicate statistically significant 
correlation / anti-correlations. 
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Percentage 


of Multi-Pulsed GRBs Showing 


Overall Pulse Behavior 
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0.54%t 
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hardness 








0.57%t 


0.59%t 
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Table 3: Multi-pulsed Short and Long GRBs have been combined together to provide a larger 
sample; least squares is applied to the properties of all pulses fitted within each multi-pulsed 
GRB to learn more about the intra-pulse distribution. The resultant gradient (slope) of each 
pulse property relative to the others indicates whether the properties within the burst are 
correlated (positive gradient) or anti-correlated (negative gradient; indicated by a f). The 
percentage of multi-pulsed GRBs exhibiting correlated behaviors is then compared to the 
results of the inter- pulse duration shown in Figure 2 and Table 2; if somewhat more than 
50% of the multi-pulsed GRBs have the same general gradient as that obtained from the 
inter-pulse distribution then the two distributions have the same general behavior. Since 
these percentages average 62%, the behaviors of the intra-pulse distribution have the same 
sense but otherwise appear to be generally weaker than the power-law forms of the inter-pulse 
distribution. 
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Table 4: Hard-to-soft (II2S) and Intensity tracki ng (IT) pulses in c ommo n with BATSE puls e 
data base. Pulse cla ssification sources are Pl= (|Peng et al.ll2009al ). P2= (|Peng et al.ll2009bl ). 
L1 = JLu et alJboioh . 
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Fig. 1. — Sample 4-channel pulse fits for BATSE triggers 1221 (upper left), 432 (upper 
right), 1443 (lower left) and 1468 (lower right). Overlaying the 64- ms data are vertical lines 
defining each Bayesian block segment, as well as the pulses finally fitted to the data. 
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Fig. 2. — Pulse correlations for Long (red) and Short (blue) GRB pulses, showing the strong 
correlation between these parameters. 
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Fig. 3. — Pulse duration vs. pulse peak lag for Long (diamond) and Short (*) GRB pulses, 
showing the strong correlation between these parameters. Only pulses with positive lags are 
plotted. 
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Fig. 4. — Pulse 256 ms peak flux vs. duration for multi-pulsed GRBs, indicating general 
intra-pulse characteristics. Each data point identifies the mean peak fiux vs. mean duration 
for a mutli-pulsed GRB; the error bars reflect the mean error. The shorter the average pulse 
duration is in a burst, the brighter the average pulse is. 
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Fig. 5. — Pulse duration vs. lag for multi-pulsed GRBs using the same notation as Figure 
4 (only bursts with positive mean pulse lags are plotted). The greater the average pulse 
duration is in a burst, the longer average lag the pulse will have. 
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Fig. 6. — Pulse hardness vs. duration for multi-pulsed GRBs, using the same notation as 
Figure 4. The harder the average pulse is for a multi-pulsed GRB, the shorter the average 
pulse will be. 
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Fig. 7. — Pulse fluence vs. duration for multi-pulsed GRBs, using the same notation as 
Figure 4. The greater average energy fluence a pulse has, the longer the average pulse will 
be. 
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Fig. 8. — Pulse 256 ms peak flux vs. duration for Long (diamond) and Short (*) GRB pulses 
for single-pulsed GRBs triggering on the 256 ms timescale. 
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Fig. 9. — Estimated limitations on pulse sampling threshold for Long (diamond) and Short 
(*) GRB pulses in terms of fluence/duration and 64 ms peak flux. Sampling is relatively 
complete above and to the right of the solid line, with the noted exception of pulses in fields 
too crowded to cleanly fit. The estimated lower fimit (below which almost no pulses have 
been fitted) is indicated by the dotted line. Approximate durations of the plotted pulses are 
given. 
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Fig. 10. — GRB pulse spectral evolution, as seen from the fits to single-pulsed BATSE 
Trigger 1883. The solid line denotes the channel 4 fit, the dotted-dashed line is channel 3 
fit, the dashed line is channel 2 fit, and the dotted line is channel 1 fit. 
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Fig. 11. — Approximat e isotropic luminosity p er pulse for Long (diamond) and Short (*) 
GRB pulses (data from iGhirlanda et al.l ( l2009l )). demonstrating that Short GRB pulses are 
not systematically less luminous than Long GRB pulses, even though Short GRBs are less 
luminous than Long GRBs. 



